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Abstract. We present and analyze quiescent UBVRI light curves of the classical symbiotic binary YY Her. We 
show that the secondary minimum, which is clearly visible only in the quiescent VRI light curves, is due to 
ellipsoidal variability of the red giant component. Our simple light curve analysis, by fitting of the Fourier cosine 
series, resulted in a self-consistent phenomenological model of YY Her, in which the periodic changes can be 
described by a combination of the ellipsoidal changes and a sinusoidal changes of the nebular continuum and line 
emission. 
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1. Introduction 
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Symbiotic stars are long-period interacting binaries made 
up of an evolved red giant and a hot companion sur- 
rounded by an ionized nebula. The hot component in 
the vast majority of systems seems to be a luminous 
(~ 1000 Lq) and hot (~ 10 5 K) white dwarf powered by 
thermonuclear burning of the material accreted from its 
, companion's wind. Depending on the accretion rate, these 
systems can be either in a steady burning configuration or 
undergo hydrogen shell flashes, which in many cases can 
last for decades and even centuries due to low mass of the 
white dwarf. The latter can be the case for non-eruptive 
symbiotics (such as RW Hya or V443 Her), in which the 
hot components maintain roughly constant luminosity and 
temperature for many decades (Mikolajewska 1997 ). In 
many systems, however, the hot components show mul- 
tiple outburst activity with timescales of a few/several 
years, which cannot be simply accounted for by the ther- 
monuclear models. Among them, Z And is one of the best 
studied (e.g. Mikolajewska & Kenyon 1996) , and references 
therein), and the symbiotic stars with multiple outburst 
activity are often referred to as Z And-type systems. 

The photometric history of YY Her sin ce 1890 has been 
studied in detail by Munari et al. ( 1997a , hereafer M97a) 
who revealed, in addition to Z And-type multiple outburst 
activity, a periodic fluctuations with P = 590 d and a vi- 



sual amplitude of 5s 0.3 m . The orbital nature of this pe- 
riodicity has been confirmed by subsequent UBV photo- 
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metric observations (Munari et al. 1997b, Tatarnikova et 
al. 2001, hereafter M97b and T01, respectively). In par- 
ticular, the photometric and spectroscopic study of YY 
Her during the major outburst in 1993 and its decline 
(Tatarnikova et al. |2000i hereafer TOO; T01) has shown 
that the minima are due to obscuration of the hot com- 
ponent and the ionized nebula. The return to quiescence 
in 1997/98 was accompanied by significant changes in the 
shape of light curves, and in particular by the appearance 
of a secondary minumum in the VRI light curves (T01; 
Hric ct al. |2001 , hereafter H01). T01 interpreted the sec- 
ondary minimum in terms of ellipsoidal changes of the red 
giant, whereas H01 argued that the secondary minimum 
is caused by an eclipse of the red giant by a circumstellar 
envelope around the hot component. 

In this paper we present new UBVRI photometric 
observations obtained during the secondary minimum in 
2001, and demonstrate that the quiescent light curves of 
YY Her can be succesfully reproduced by a combination of 
ellipsoidal changes of the cool giant and sinusoidal changes 
of the nebular continuum and line emission. 



2. Observations 

UBVR'I 1 photometry was obtained with the 0.6-m and 
1.25-m telescopes at the Crimean Station of the Sternberg 
Astronomical Institute (SAI). For UBV bands (Table y 
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Table 1. UBV photometry of YY Her 



Table 2. BVR'I' photometry of YY Her 
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the 2- channel photometer with an S-20 photomulti- 
plier was used, wheras a CCD matrix was the detec- 
tor for BVR'I' bands (Table ^; see below for the defi- 
nition of R'l'). The uncertainty on individual measure- 
ments does not exceed 0.03 mag. The UBV bands corre- 
spond to Johnson's broad-band system. The instrumental 
ri bands have the effective wavelengths, A r ss 7000 A and 
A; w 9000 A, and the band widths AA r w 2400 A and 
AAj « 3000 A, respectively, and they are different from 
any commonly used photometric system. The ri observa- 
tions were reduced to Johnson's RI magnitudes using the 
standard stars from M67 (Mendoza, 1967). In the case of 
stars with a combination spectra such as YY Her, a trans- 
formation of the instrumental magnitudes to the standard 
system may be not as precise as in the case of normal stars. 
For this reason, we refer to these magnitudes as R'l'. HD 
168957 (U = 6 m 35, B = 6 m 91, V = 7 m 01), and the star 
designated as G (V = 13. m 08, R' = 12. m 44, V = ll. m 94) 
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Fig. 1. The UBVRI photometry of YY Her in 1993-2001. 
Closed symbols represent our UBV R'l' data (this study; 
M97a; TOO; T01), and open symbols the B C V C R C I C data 
of H01) transformed to our system (see text). 



on the finding chart in M97a) were adopted as the com- 
parison stars. 

Our data have been combined with published UBV 
photometry (T01, and references therein) as well as 
B C V C R C I C one performed in the modified Johnson-Kron- 
Cousins system by H01). The resulting light curves are 
shown in Fig. [j]. The B C V C R C I C magnitudes in Fig. [I] 
were transformed to our system assuming B = B c + 0.3, 
V = V c + 0.2, R' = R c ~ 0.2, and I' = I c — 0.52, respec- 
tively. 



3. Analysis 

The most prominent features of the UBV light curves of 
YY Her (Fig. |l|) are deep eclipses and the large outburst 
that started in 1993. The R C I C /R'I' light curves cover the 
late decline and return to quiecence, and they show the 
characteristic double hump. The secondary minimum is 
also clearly visible in the quiescent V light curve as well 
as the most recent B data. 

In our analysis, we will focus on the quiescent light 
curves. Fig. shows the UBVRI magnitudes folded 
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over the last three orbital periods. We have adopted the 
ephemeris, 



Min = JD 2 450 686.2(±2.1) + 589.5(±0.3) x E, 



(1) 



derived from 13 individual minima estimates (by a 
parabolic fit) in published VRI light curves (M97a; M97b; 
TOO; TOf; HOI). The salient characteristic of the quies- 
cent light curves, namely: (i) nearly sinusoidal shape of 
the U light curves; (ii) a flat toped B light curve; and (iii) 
the two minima of almost comparable depth in I 1 light, 
combined with the fact that the quiescent, ultraviolet + 
optical, spectra of YY Her can be satisfactorily matched 
with a three-component model composed of a cool giant 
dominating the red spectral range, a gaseous nebula dom- 
inating the near UV range, and the hot component pre- 
dominant in the short UV flux (TOO), led us to a simple 
phenomenological model in which the light changes are 
described by a combination of variable nebular emission 
and ellipsoidal changes of the red giant. 

Each light curve has been converted to fluxes, F\ = 
10 _0 4mA , and then independently fitted by the Fourier 
cosine series 



F\ = Fx,o - Ai,a cos 4> - A g ^x cos 2</>. 



(2) 



The results of our fits to the UBVR'I' (our data) and 
R C I C (HOI data) light curves are summarized in Table ||, 
columns 2-4. In the case of U data, the best fit is ob- 
tained with the cos<ji> term only. This indicates that the 
contribution of the giant to the U light is negligible, 
and that all changes can be attributted to the nebu- 
lar emission. Assuming that the amplitude of the neb- 
ular emission changes is the same at all wavelenghts, 

{Fix + A*M*Zx ~ A n,x) « (Fv.o + Av)/(F Vj0 - A v ) = 
C ~ 4.5, we can separate the contributions of the nebular 
and the giant components to the constant F\£. The re- 
sulting charateristics of the nebular source (columns 5-7) 
and the red giant (columns 8-10) are also given in Table |[ 



The maximum, Fx 



0-5)=F°+A ni x, fluxes 



of the nebular component have been calculated using the 
zeroth magnitude star fluxes from Straizys ( 1992 ). The 
UBVR'I' and R C I C magnitudes corresponding to these 
maximum fluxes have been also calculated. We note here 
that the broad-band magnitudes, and the resulting fluxes 
may be significantly afftected by emission lines. Using the 
quiescent optical/red spectra of YY Her (TOO; T01), we 
have calculated the contribution of the emission lines to 
the observed BVR' magnitudes as Am £ 0.15/0.02/0.08 
mag, respectively (this contribution can be, in fact, even 
higher because we have not included faint lines). Using 
these estimates, we have found that the nebular line emis- 
sion is responsible for at least 30%(£? band), 7%(V) and 
50 %(-/?') of the total nebular emission fluxes, and the cor- 
responding magnitudes, given in the 5th and 6th columns 
of Table Finally, the maximum emission measure of the 
nebular continuum (7th column of Table ||) has been es- 
timated from the maximum, Fx max, fluxes of the nebular 
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Fig. 2. Phased quiescent light curves of YY Her. The same 
symbols as in Fig. [l| are used. Our best fit models (Tabic ^j) 
are plotted as solid lines (our data) and dashed lines (H01 
data), respectively. 



component corrected for the contribution from emission 
lines and the interstellar reddening £b-v = 0-2 (M97b), 



EM = n 2 c V/4nd 2 = F A , max /e* f+ff . 



(3) 



c bf+ff 
A ) 



mated using our UBVR'I' filter transmission functions 
and bf+ff continuum emission coefficients (case B recom- 



bination) from Pottash ( 1984 ). The EM values derived 
for each band agree with each other to within a fac- 
tor of 2, and they cluster around EM ~ 2 x 10 14 and 
- 3 x 10 14 cm- 5 for T c = 10 4 and 2 x 10 4 K, respectively. 

The nebular emission measure derived from our light 
curve analysis requires a Lyman continuum luminosity of 
~ 6-5 x 10 45 (d/lkpc) 2 phots" 1 for T c = 10 4 -2 x 10 4 K, 
which corresponds to a hot component with Th ~ 10° K 
and Lh ~ 80(d/l kpc) 2 LQ. This can be compared with 



4 



J. Mikolajewska et al.: Tidally distorted giant in YY Herculis 



independent estimates for the quiescent L^. In particu- 
lar, the bolometric fluxes for the hot component given by 
M97b and TOO correspond to L h ~ (30-60) (d/1 kpc) 2 L . 

The UBVR'I 1 and R C I C magnitudes of the cool com- 
ponent have been calculated using to,q = — 2.5 log-Fg A , 



-2.51og(F" +A, A ),andm 



-2.51og(F° 



^.g,A), where F® A is the cool component contribution to 



the constant F\ t o- The colours derived from these magni- 
tudes, B — V = 1.75/1.77/1.73, V — R' = 1.81/1.73/1.92 
and V — I' = 3.45/3.33/3.61 representing the average, 
maximum and minimum values, respectively, are in ex- 
cellent agreement with those expected for a moderately 
reddened M3-M4 giant. Moreover, they indicate some- 
what higher temperature (earlier spectral type) at maxima 
{4> = 0.25 and 0.75) than that at minima (</> = and 0.5) 
as expected for a tidally distorted giant. 

Assuming the Roche-lobe filling factor i? g /i?L ~ 1, 
the gravity-darkening exponent, a = 0.32 (convective en- 

o 

velope) , and the system inclination i ~ 70-90 , the ampli- 
tude of the cool component changes in V is consistent with 
a mass ratio q — M g /M^ ~ 2-3. This estimate is based 
on synthetic light curves calculated using the Wilson- 
Dcvinncy (WD) code (see Belczyhski & Mikolajewska 



199S for details). The amplitudes in R, V bands, however, 



require a much higher value of a iS 1. Similar effect was 
found in T CrB, and it is most likely due to the blackbody 
approximation for wavelength dependence in the WD code 
used to calculat e the synthetic light curves (Belczyhski 
& Mikolajewska 1998 ). This assumption is probably not 
valid in the case of M-type stars with strong contribution 
of TiO bands in the optical and red p art of the spectrum. 
Recently, Orosz & Hauschildt (|2000|) have demonstrated 
that, in general, the light curves for tidally distorted cool 
giants computed using the NEXTGEN model atmospheres 
for cool giants have deeper minima in the optical range 
then their blackbody counterparts. 

Finally, we note that the amplitude of ellipsoidal vari- 
ability is a strong function of t he Ro che-lobe filling factor, 
A to cx {R g /R L ) 3 (e.g. Hall |l990| ). The lower limit for 
the cool component changes set by the uncertainty of our 
fits (Table |) then requires R g /R L ^ 0.9 for q Z 1 and 

i > 70°. 

The presence of a Roche-lobe filling giant provides an 
opportunity to estimate the distance d to YY Her. The 
Roche lobe radius is Rl/o. ~ 0.4 4±0.0 5 for the mass ratio 
of q = M g /Mh ~ 1-3 (Paczyhski 1971 ). Our orbital period 
requires the binary separation, a w 1.85 ± 0.15 a. u. for a 
reasonable total m ass of the system, M g + Mh ~ 2-3 Mq 
(e.g. Mikolajewska |1997[), so R L g 0.8±0.2a.u. Using the 



Barnes-Evans relation (Calm |l980 ) 



F K = 4.2207 - 0.1K - 0.5 log s, 



(4) 



with the K surface brightness, F K = 3.84 (M4III), and 
Kq = 7.9, we estimate the giant's angular diameter, s — 
0.152 mas, and d = 10 ± 3kpc. Such a distance implies a 
very high ~ 3 kpc height above galactic plane and locates 
YY Her in galactic halo. 



4. Concluding remarks 

Our very simplistic light curve analysis, namely indepen- 
dent fitting of the Fourier cosine series to the quiescent 
UBVR'I' light curves, resulted in a self-consistent phe- 
nomenological model of YY Her, in which the periodic 
changes can be described by variations of the nebular 
emission along the orbit (the cos <fi term) combined with 
the ellipsoidal variability of the red giant (the cos 2<fi term). 

YY Her is not the only symbiotic system with the 
secondary minima. Similar light curve behaviour is ob- 
served in CI C yg, w hose quiescent UBVRI light curves 
(Mikolajewska 2001) are almost identical with the light 
curves shown in Fig. ^. The secondary minima are also 
present in quiescent near infrared light curves of BF Cyg 
(Mikolajewska et al. |M§). 

We also think that the alternative interpretation of the 
secondary minimum in terms of the eclipse of the cool gi- 
ant by the circumstellar envelope of the hot component 
(H01) is implausible. First, it cannot be the ionized nebu- 
lar envelope because to cause eclipses it must be optically 
thick in the optical/red continuum whereas the presence of 
a prominent Balmer jump indicates that it is not. Second, 
to account for the amplitudes of the secondary minimum, 
it must screen ~ 20% of the red giant's projected sur- 
face. Such a big and hot screen, however, should dominate 
the quiescent optical continuum, which is not observed. 
Finally, we can also exclude eclipses of the cool giant by 
an accretion disc. A geometrically and optically thick disc 
seen nearly edge-on could mimic a cool source, and screen 
significant part of the giant's surface. Such a disc, however, 
would also completely hide the hot white dwarf wheras the 
quiescent shortwavelength IUE spectrum clearly shows a 
hot continuum (M97b). Similarly, a hot rising continuum 
is present in quiescent shortwavelength IUE spectra of CI 
Cyg and BF Cyg. 

One of fundamental questions in relation of symbiotic 
binaries is the mechanism that powers the Z And-type 
multiple outburst activity observed in many classical sym- 
biotic systems. The quiescent spectra of these systems can 
mostly be fitted by a hot stellar source (most likely a white 
dwarf powered by thermonuclear burning of the material 
accreted from its companion) and its ionizing effect on 
a nebul a (e.g . Miirset et al. |l99l| ; M97b, Mikolajewska & 
Kenyon 1996). Their outburst activity, however, with time 
scales of a few/several years cannot be simply accounted 
by the thermonuclear models. A promising interpretation 
of this activity involves changes in mass transfer and/or 
accretion disc instabilities. Detection of an ellipsoidal hot 
continuum source during outbursts of CI Cyg, AX Per, YY 
Her, AS 338 and other Z And-type systems suggests the 
presence of an optically thick accretion disc (Mikolajewska 
& Kenyon |1992|; Esipov et al. 120001 TOO; Mikolajewska et 



al. p002[) strongly supports this interpretation. 

Related to this problem is the process of mass transfer 
- Roche lobe overflow or stellar wind - and the possibility 
of an accretion disc formation. The presence of secondary 
minima in the light curves of YY Her, apparently due to 
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Filter Best fit parameters Hot nebular source Cool giant 
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n, A 


A™ 


F ( 2 ) 

A, max 


m max [mag] 


EM (3) 


mo [mag] 


m max [mag] 


m min [mag] 
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2.65 ±0.05 


1.68 ± 0.06 





17.8 ±0.4 


13.41 ±0.03 


1.5/2.7 








B 


1.77 ±0.02 


0.54 ± 0.03 


0.17 ±0.03 


8.6 ±0.3 


14.64 ± 0.06 


3.0/2.4 


15.21 ±0.06 


15.03 ±0.08 


15.43 ±0.11 


V 


5.64 ± 0.05 


0.96 ± 0.06 


0.84 ±0.07 


9.0 ±0.4 


14.02 ±0.06 


2.8/3.2 


13.46 ± 0.04 


13.26 ±0.05 


13.70 ±0.07 




26.4 ±0.4 


2.9 ±0.4 


2.5 ±0.4 


12.4 ±0.9 


12.82 ±0.14 


2.0/2.9 


11.65 ±0.03 


11.53 ±0.04 


11.78 ±0.05 




26.5 ±0.2 


3.4 ±0.3 


1.7 ±0.3 


19.7 ±0.9 


12.64 ±0.09 


3.2/4.5 


12.00 ±0.03 


11.89 ±0.04 


12.13 ±0.05 


J'(4) 


101.9 ±0.9 


1.5 ±0.7 


7.2 ± 1.1 


3.0 ±0.9 


13.5 ±0.5 


1.6/2.0 


10.01 ± 0.02 


9.93 ± 0.03 


10.09 ± 0.04 


r(5) 


63.5 ±0.6 


1.6 ±0.8 


5.4 ±0.8 


4.9 ± 1.5 


13.5 ±0.5 


2.5/3.0 


10.54 ± 0.03 


10.44 ± 0.04 


10.64 ± 0.05 



(1) in units of 10 6 ; (2) in units of 10 15 ergs 1 cm 2 A 1 ; (3) maximum emission measure in units of 10 14 cm 5 , for 
T c = 10 4 /2 x 10 4 K; (4) our R'/I' data; (5) R c /L from H01). 



ellipsoidal changes of the red giant, as well as in a few 
other similar systems provides important clues. Although 
it is very premature to claim that all symbiotic systems 
with the Z And-type activity do have tidally distorted 
giants, and - at least during active phase - accretion discs, 
whereas the non-eruptive systems (such as RW Hya) do 
not, the former is certainly the case for YY Her, CI Cyg 
and a few other active systems. 

Systematic searches for ellipsoidal variations in both 
active and non-eruptive symbiotic stars are necessary to 
address the problem, and confirm or exclude tidally dis- 
torted giants. The observations must be, however, carried 
in the red and near-IR range where the cool giant dom- 
inates the continuum light. We also note that any opti- 
cal/red light curve analysis must be supported by spec- 
troscopic information about contamination of the broad- 
band photometry by the nebular component. 
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